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Abstract 

As an alternative to dark energy it has been suggested that we may be at the center of an inhomo- 
geneous isotropic universe described by a Lemaitre-Tolman-Bondi (LTB) solution of Einstein's field 
equations. In order to test this hypothesis we calculate the general analytical formula to fifth order 
for the redshift spherical shell mass. Using the same analytical method we write the metric in the 
light-cone by introducing a gauge invariant quantity G(z) which together with the luminosity distance 
Dl(z) completely determine the light-cone geometry of a LTB model. 
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I. INTRODUCTION 



High redshift luminosity distance measurements 
measurement 



24 



32 



33 



34 



35 



37 



36j and the WMAP 



39] of cosmic microwave background (CMB) interpreted in the context of 



standard FLRW cosmological models have strongly disfavored a matter dominated universe, and 
strongly supported a dominant dark energy component, giving rise to a positive cosmological 
acceleration, which we will denote by a FLRW (not to be confused with the scale factor a). As 
an alternative to dark energy, it has been proposed that we may be at the center of 

an inhomogeneous isotropic universe described by a Lemaitre-Tolman-Bondi (LTB) solution 
of Einstein's field equations, where spatial averaging over one expanding and one contracting 
region is producing a positive averaged acceleration ap. Another more general approach to 
map luminosity distance as a function of redshift Dl(z) to LTB models has been recently 



proposed [27|, 



281 ] . showing that an inversion method can be applied successfully to reproduce 



the observed Dl(z). 

The main point is that the luminosity distance is in general sensitive to the geometry of 
the space through which photons are propagating along null geodesies, and therefore arranging 
appropriately the geometry of a given cosmological model it is possible to reproduce a given 
Dl(z). For FLRW models this corresponds to the determination of Q\ and Q m and for LTB 
models it allows to determine the functions E(r), M(r),tb(r). 

Another observable which could be used to constraint LTB models is the redshift spherical 
shell mass mn(z), which we calculate for a central observer up to the filth order in the red-shift. 
Using the same set of geodesic equations derived to obtain such central expansion we also write 
the light-cone metric in red-shift space. 



II. LEMAITRE-TOLMAN-BONDI (LTB) SOLUTION 



iflB 



Lemaitre-Tolman-Bondi solution can be written as 

ds* = -dP + + R 2 cM 2 , (1) 

1 + 2 hi 

where R is a function of the time coordinate t and the radial coordinate r, R = R(t,r), E is 
an arbitrary function of r, E — E(r) and R, r = dR/dr. 
Einstein's equations give 

fR\ 2 _ 2E(r) 2M(r) 

\R ~ R 2 R 3 ' { > 



p{t,r) 



Myr 

R 2 R, T 



(3) 



with M = M(r) being an arbitrary function of r and the dot denoting the partial derivative 
with respect to t, R = dR(t,r)/dt. The solution of Eq. ([2]) can be expressed parametrically in 
terms of a time variable r = /' dt'/R(t', r) as 

M(r) 



Y r,r) 



t(r,r) 



-2E(r) 

M(r) 
-2E(r) 



1 - cos ( y—2E(r)r 
1 



-2E(r) 



sm 



-2E(r)r 



tb(r) 



(4) 
(5) 



where Y has been introduced to make clear the distinction between the two functions R(t, r) 
and Y(t, r) which are trivially related by 



R(t(r,r)) = Y(r,r) 



(6) 



and tfy(r) is another arbitrary function of r, called the bang function, which corresponds to the 
fact that big-bang/ crunches can happen at different times. This inhomogeneity of the location 



of the singularities is one of the origins of the possible causal separation 
observer and the spatially averaged region for models with positive ar>. 
We introduce the variables 

2E{r) 6M(r) 



29] between the central 



A(t,r) = ^^, k(r) 



Po{r) 



(7) 



so that Eq. (OQ) and the Einstein equations (T5]) and (131) are written in a form similar to those 
for FLRW models, 



rfs 2 = _ dt 2 + A 2 



1 + 



A, r T 



A 



dr 2 



1 — k(r)r' 



+ r 2 dn 2 2 



(8) 



The solution of Eqs. (j4j) and 

a(r],r) = 



(3 

p(t,r) 



H r ) , Po(r) 



A 2 
(p r 3 ) 



3A 3 



6A 2 r 2 (Ar), 
]) can now be written as 
Po{r) 



6k(r) 

Po(r) 
6k(r) 



1 — cos \\Jk{r) rj 



sm 



k(r) 



k{r) Tj 



+ t b (r) 



(9) 
(10) 



11) 



12) 
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where r] = rr = /* dt' / 'A(f , r) and A(t(7p r), r) = a(rp r). 

In the rest of paper we will use this last set of equations . Furthermore, without loss of 
generality, we may set the function po(r) to be a constant, po(r) — Po — constant, corresponding 
to the choice of coordinates in which M(r) oc r 3 , and we will call this, following JS], the FLRW 
gauge. 

We need three functions to define a LTB solution, but because of the invariance under general 
coordinate transformations, only two of them are really independent. This implies that two 
observables are in principle sufficient to solve the inversion problem of mapping observations to 
a specific LTB model, for example the luminosity distance Dl{%) an d the redshift spherical shell 



mass m(z)n(z) = mn(z). As observed by 17|, there as been sometime some confusion about 
the general type of LTB models which could be used to explain cosmological observations, so 
it is important to stress that without restricting the attention on models with homogeneous 
big bang , £&(r) = 0, a void is not necessary to explain both Dl(z) and m(z)n(z) = mn(z) 
without cosmological constant. We call mn(z) redshift spherical shell mass, since this quantity 
it is not the galaxy number counts as it is called in the interesting paper [17|, but the product 
of the source number density n(z) times the source mass function m(z), and it has dimension 
of energy. 



This should not be confused with the redshift spherical shell energy Erss introduced in 41] , 
which is a quantity obtained by integrating mn(z) over varying redshift intervals AZ(z) 



z+AZ(z) 

E RSS (z)= J 4nmn(z')dz' (13) 

z 

t{z)-t(z + AZ(z)) = At (14) 

corresponding to the same constant time interval At, which, if chosen to be sufficiently 
smaller than the time scale of astrophysical evolution of the source, should eliminate the effect 
of the source evolution on mn(z). 

III. GEODESIC EQUATIONS 

We will adopt the same method developed in [ijj] to find the null geodesic equation in the 
coordinates (i],t), but here instead of integrating numerically the differential equations we will 
find a local expansion of the solution around z = corresponding to the point (to, 0) = (rjo, t), 
where to = t(voi r )- We will also provide more details about the geodesic equation derivation 



which were presented in |14|] in a rather concise way. We will indeed slightly change notation to 
emphasize the fully analytical r.h.s. of the equations obtained in terms of (77, t), on the contrary 
of previous versions of the light geodesic equations which require some numerical calculation of 
R{t,r) from the Einstein's equation (j2J). 

For this reason this formulation is particularly suitable for the derivation of analytical results. 

The luminosity distance for a central observer in a LTB space as a function of the redshift 
is expressed as 

D L (z) = (1 + zfR (t(z),r(z)) = (1 + z) 2 r(z)a {r){z) : r{z)) , (15) 

where (t(z),r(zfj or ((t](z), r(z) J is the solution of the radial geodesic equation as a function 
of the redshift. 

The past-directed radial null geodesic is given by 

dT(r) -R r (t,r) 

— — = /(T(r),r); f{t,r) = —= 16 
dr <'l + 2E(r) 



where T(r) is the time coordinate along the null radial geodesic as a function of the the coor- 
dinate r. 

From the implicit solution, we can write 

T(r) =t(U(r),r) (17) 

dr(r) = mdu<j)_ + m 

dr dr\ dr dr 

where U(r) is the 77 coordinate along the null radial geodesic as a function of the the coor- 
dinate r. 

Since it is easier to write down the geodesic equation in the coordinate (t, r) we will start 
from there 
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dr y/l + 2E(r(z)) 

dz (l + z)R'[r{z),t{z)\ [ ' 

dt _ R![r{z),t{r)\ 
~d~z ~~ ~ (l + z)R'[r(z),t(z)]' 

(20) 

where the ' denotes the derivative respect to r and the dot ' the derivative respect to t. These 
equations are derived from the definition of redshift and by following the evolution of a short 
time interval along the null geodesic T(r). 



The problem is that there is no exact analytical solution for R(t,r), so the r.h.s. of this 
equations cannot be evaluated analytically but requires to find a numerical solution for R first 



26] , and then to integrate numerically the differential equation, which is a quite inconvenient 
and difficult numerical procedure. 

Alternatively a local expansion for R(t,r) around (to;0) , corresponding to the central ob- 
server, could be derived and used in eq. (l20l) . but being an expansion will loose accuracy as the 
redshift increases. 

For this reason it is useful for many numerical and analytical applications to write the 
geodesic equations for the coordinates (r],r), 

dr] d r t(r],r) — F(r],r) 



dz (1 + z)d ri F(r], r) 

dr G-(?7, r ) 

dz (1 + z)d v F(r], r) 



p(rj,r), (21) 
qfar), (22) 



F(r),r) = , 1 [d r (a(rj,r)r) + d v (a(r),r)r)d r r)] , (23) 

J 1 — k(r)r 2 

where rj = U(r(z)) and F(r),r) = f (t(r] , r) , r) . It is important to observe that the functions 
p, q, F have an explicit analytical form which can be obtained from a(r), r) and t(r), r) as shown 
below. 

The derivation of the implicit solution a(r], r) is based on the use of the conformal time 
variable rj, which by construction satisfies the relation, 

drj(t, r) 



dt 

This means 



a' 1 . (24) 



t(r},r) = t b (r)+ / a(rj ,r)dr) , (25) 
Jo 

dt = a(r}, r)dr] + (£ ^|^W + t b (r)j dr , (26) 

In order to use the analytical solution we need to find an analytical expression for F and F^. 
This can always be done by using 



d r t{7],r) 



3 sin (v\fk(r)) — rj ( 2 + cos {rj\Jk{r) ) \J k{r) 



I2k(rf/ 2 L" V v '" v ' V '' V" ' V v " v ' ) v " v ' ) J + ' ^ 
d r r, = -a(r ll r)' 1 d r t (29) 

In this way the coefficients of equations ( I2T1) and (J22l) are fully analytical, which is a significant 
improvement over previous approaches. 



IV. LIGHT-CONE METRIC 



Using this equations we can write the metric on the light-cone in redshift space by using 
dr = (dr/dz)dz and 



dt , / dt dr] dt dr\ , . , 

dt = —dz = + — — dz 30 
dz \ or] dz dr dz J 

p(z) = p(r](z),r(z)) (31) 

q(z) = q(r](z),r(z)) (32) 

R(z) = R(r)(z),r{z)) = r(z)a(rj(z),r(z)) (33) 



to finally get 



ds\ c = G{z)dz 2 + R(z) 2 dQ 2 (34) 
G(z) = q(z) 2 F(r](z),r(z)) 2 -Q(z) 2 (35) 
Q{z) = p{z)d v t{r]{z),r{z)) + q{z)d r t{r){z),r{z)) (36) 

from which we can find directly the radial null geodesic equation for rj(z) by imposing 
G(z) =0. This metric only describes space-time in the light-cone, and it is not valid outside, 
since it is based on the null geodesic congruence given by the radial null geodesies tangent 
vector field. Since this is the only part of the full LTB space observationally connected to 
a central observer, it is sufficient to determine models which can explain observational data, 
which by definition have to be inside the light cone of the central observer. 

In principle, given G(z) and R{z) = Da{z) the light-cone geometry of a LTB is completely 
defined since this are gauge invariant quantities, but it is not clear what observable to associate 
to G(z), so the inversion problem still requires to solve the geodesic equations to construct 
mn{z) and D L (z) = (1 + z) 2 Da{z) . 

It should be mentioned that our approach is general while Q derived an analytical version of 
the geodesic equation in the light-cone gauge, t —t = —r, under the assumption E{r) = 0, just 
taking into account the inhomogeneity coming from the bang function £&(r). Another related 
result was obtained by , working in the same gauge, where an expression was derived to third 
order in red-shift for r(z) and Dl(z). 

Our results are equivalent after performing the appropriate gauge transformations obtained 
from the condition 
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M(f) = M(r) . (37) 

Such transformations involve the relation between the coefficients of the expansion of 
M(f), E(r) and our k(r) = —2E(r)/r 2 , which are rather complicated already at third or- 
der in red-shift. For this reason we will report them separately in a future work in which the 
general relation between the light-cone gauge and the FLRW gauge is analyzed in detail. 



V. CALCULATING mn{z) 

Expanding the r.h.s. of the geodesies equation we can easily integrate the corresponding 
polynomial in q(z),p(z), to get r(z) and r](z). It can be easily shown that in order to obtain 
Dl(z) to the fourth order and mn(z) the fifth we need to expand r(z) to the fourth order and 
T](z) to the third. 

Since we only need r(z) to calculate mn(z) we only give that expansion, but it is understood 
that r](z) has to be computed as well to obtain it. 

In order to have an solution which is analytical everywhere we will should use the following 
expressions for k(r) and t&(r): 



k(r) = k + k 2 r 2 + hr 4 (38) 
t b {r) = t b + t b 2 r 2 (39) 

which are based on the fact that taking only even powers the functions are analytical ev- 
erywhere, including the center. We will nevertheless keep in the following formulas, in order 
to provide the general results and show where the differences arise respect to the previous 
calculations 8] to lower order in redshift. 

After re-expressing the results in terms of H and go we get 



r \Z + r 2 z 2 + r 3 z 3 + r^z 4 
1 

go + 1 
2H 

(2q - 1) (# 4 (1 - 2q ) 2 f-^t + q 2 + ^ _ 5^ + + Qk 2v ^=Tq 2 arccos (^) 

2H$(2q - 1)3 
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r 4 



# 4 (1 - 2g ) 2 (5g 3 - q 2 + 4) - 2ff 3 (l - 2g ) 2 g (10g + l)t b 2 - 3A; 2 (10g 2 + 5q - 2) 

8# 5 (2g - l) 2 

8# 2 g (l - 2q )H b + 6k 2 q 2 ^2q^l(10q + 1) arccos 



8H$(2q - I) 3 



(40) 



where 



H. = 



qo = - 



g(t ,0) 
a(t ,0) 

a(f , 0)d(t o , 0) 



(41) 
(42) 



a(to,0) 2 

The derivative respect to t is denoted with a dot, and is calculated using the analytical 
solution a(i],r) and the derivative respect to r] is obtained from d t a — a — d^aa^ 1 . 
From the definition of mn(z) and the equation for the energy density we can write 



4:7imn(z)dz = pd 3 V 
from which by using dr = (dr/dz)dz we get 



AtxM 1 



-.dr 



mn(z) 



M'Mz) 



dr(z) 



\Jl — k{r)r- 



p r\z) 



(43) 



dr(z) 



(44) 



yjl - k(r(z))r(z) 2 dz 2yjl - k(r(z))r(z) 2 dz 
where in the last equation we have used the FLRW gauge condition M(r) = p r 3 /6, which 
allows to calculate mn{z) directly from r(z). 
We finally get: 



mn{z) 



3q z 2 6q (q +l)z 3 



3g 



4/J 5 (-l + 2g )(5/2) 
+20Fq( 1 - 2go) 2 goV / 2?o - U b 2 + Q0k 2 q 2 arccos 



2q - 1 (# 4 (1 - 2g ) 2 (l5g 2 + 14g + 13) + 10fc 2 (l - 5q j) 

1 N 



Wo 



+ 



y/2q - 1 (#o(l - 2<?o) 2 (28g 3 + 24g 2 + 21g + 19) - 3A; 2 (50g 2 + 31g - 10 



+^q ~ 1 (24F 2 (1 - 2q ) 2 q t b 3 - 6# 3 (1 - 2q ) 2 q (Uq + 5)4) 

1 N 



+18/c 2 (14g + 5)g arccos 



(45) 



As it can be seen the effects of inhomogeneities show only from the fourth order, and there 
is no dependency on fc 4 up to fifth order. 
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A similar expansion can be obtained for the luminosity distance, and then by combing the 
two results together we can give a general solution to the inversion problem at low-redshift. 

Finally it should be mentioned again that in order to avoid a central singularity we should 
set 4 = 0. 

VI. RELATION OF mn{z) TO OBSERVATIONS 

The redshift spherical shell mass mn(z) is the product of the number of sources n(z) times 
their mass m(z), so it can be related to observations by using the following trivial relation with 
the total rest mass M T (z) within a given redshift, i.e. contained in a sphere of comoving radius 



The quantity Mt(z) is obtained from observations by simply adding the mass of all the 
sources with redshift equal or less than z. 

Alternatively if we are interested in a more direct relation of mn{z) to observations, without 
having to use M T (z), E RSS (z) is the quantity which should be naturally considered, since the 
uncertainty in the redshift determination would always imply the necessity of same integration 
in redshift space of mn(z). 

The relation between mn(z) and Erss( z ) given in eq. (13,14) is particularly useful for ob- 
servational purposes, since an appropriate choice of the times scale AT can significantly avoid 
unwanted astrophysical evolution effects on the sources number counts n(z), which is exactly 
the reason why E RS s(z) it is defined in that way. 

VII. CONCLUSION 

We have derived a set of differential equations for the radial null geodesies in LTB space- 
time without a cosmological constant and applied them to compute the redshift spherical shell 
mass to fifth order in the red-shift, which could be used to test local inhomogeneities to a 
higher level of accuracy. We have also used the geodesic equations to write the light-con metric 





(46) 



(47) 
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in the red-shift space, and show clearly the relation between LTB models and cosmological 
observables. 

In the future it will be interesting to extended our analysis to other observables such as H(z) 
and w(z) and to use the equations we derived to provide a new method to solve the inversion 
problem of mapping the observed luminosity distance Dl(z) and the redshift spherical shell 
mass mn(z) to LTB models. We can in fact write the differential equation d z ( ) = 



d z (r(z)a(rj(z) , r(z)) , where the Dl{z) is the observed luminosity distance and the r.h.s. can be 
expressed analytically using the geodesic equations we obtained. This equation together with 
the one we derived for mn(z) should allow to solve for k(r(z)) = k(z) and tb(r(z)) = tb(z), with 
analytical results at low red-shift and numerical at higher red-shift. A detailed analysis of such 
approaches will presented in a future paper. 

Anther possible application of our results would be to give an analytical approximation for 
Erss{ z ) |4jJ , the redshift spherical energy, a quantity constructed by integrating mn(z) over 
varying redshift intervals AZ(z) corresponding to a constant time interval At, which should be 
the characteristic time scale over which the astrophysical evolution of the astrophysical object 
counted can be neglected. 

The formula we have derived for mn(z) could be used to test local inhomogeneities in red- 
shift space in a self-consistent way, and it will be the subject of a future paper the corresponding 
analysis of experimental data from galaxy surveys. 
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